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SECTION  I 
INTRODUCTION 

1 . %  BACKGROUND  ^  7 

y 

Fiber-reinforced  composite  materials  are  being  introduced  into  pri¬ 
mary  structures  of  modern  jet  fighter  aircraft,  e.g.,  vertical  and 
horizontal  stabilizers  on  the  F-15  and  F-16  fighters.  These  materials 
are  made  up  of  many  layers;  each  layer  consists  of  many  rows  of  parallel 
fibers  bound  together  by  a  matrix  material  of  different  composition  than 
the  fibers.  Failure  often  occurs  in  these  composites  via  delamination 
which  initiates  at  gas  bubbles  (porosity)  remaining  in  the  matrix  after 
the  manufacturing  process.  For  this  and  other  reasons,  non-destructive 
inspection  techniques  capable  of  determining  the  presence  of  the  porosity 
are  being  sought.  Current  efforts  include  inspection  with  elastic  waves. 
Items  of  interest  include  comparison  of  the  response  from  the  fibers  to 
that  from  porosity  and  the  degree  to  which  the  fibers  may  appear  trans¬ 
parent  to  the  waves..- 

For  this  study,  the  fibers  are  modeled  as  cylinders  and  the  porosity 
as  isolated  spheres.  Many  studies  have  been  done  on  the  scattering  from 
cylinders  and  isolated  spheres  (References  1-3),  linear  arrays  of 
parallel  cylinders  (Reference  4),  and  two-dimensional  arrays  of  parallel 
cylinders  (References  5,  6).  These  studies,  however,  all  provide 
frequency-domain  solutions.  Current  inspection  techniques  utilize 
broadband  time-domain  pulse  techniques;  thus,  a  time-domain  solution  for 
the  scattering  might  be  more  readily  applied  to  the  inspection  problem. 


2.  PROBLEM  AND  SCOPE 


The  problem  investigated  in  this  study  is  the  solution  for  the  dila¬ 
tation  wave  backscatter  from  cylindrical  and  spherical  inclusions  result¬ 
ing  from  incident  dilatation  pulses  which  are  representative  of  the 
pulses  produced  by  piezoelectric  transducers.  Obtaining  and  exploring 
the  differences  between  backscatter  from  spheres  and  cylinders  and 
comparison  with  an  experimental  result  were  the  major  objectives  of  this 
project.  A  secondary  objective  was  a  study  of  the  transparency  of  common 
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composite-reinforcing  fibers  to  elastic  waves  in  order  to  assess  the 
feasibility  of  non-destructive  inspection  with  backscattered  dilatation 
waves . 

The  analysis  is  limited  to  backscatter  from  single  cylindrical  or 
spherical  inclusions  of  arbitrary  homogeneous  anisotropic  elastic  compo¬ 
sition  embedded  in  a  homogeneous  isotropic  host.  The  solution  for  the 
cylinder  is  obtained  for  normal  incidence  of  the  dilatation  wave  with 
respect  to'  the  cylinder's  axis  of  symmetry.  All  media  are  considered 
to  obey  the  laws  of  linear  elasticity  and  are  assumed  to  be  non-attenua- 
tive. 

3.  APPROACH  AND  PRESENTATION 

A  review  of  Lee's  development  (Reference  7)  of  a  time-domain  integral 
equation  for  the  scattered  field  is  given  in  Section  II,  followed  by  a 
statement  of  the  time-domain  first  Born  approximation  and  the  form  of 
the  solution  for  the  scattered  pulse  resulting  from  this  approximation. 

A  "transparency  condition"  is  then  obtained  for  a  general  anisotropic 
homogeneous  inclusion  in  a  homogeneous  isotropic  host.  In  Section  III, 
the  solution  for  backscatter  from  cylinders  and  spheres  is  worked  out 
in  detail.  Power  series,  Laplace  transform  methods,  and  recursion 
relations  are  used  to  obtain  the  solutions.  Some  specific  results  in 
the  form  of  graphs  of  backscattered  responses  are  presented  and  discussed 
in  Section  IV,  where  comparison  to  an  experimental  result  is  also 
presented.  The  transparency  of  typical  fibers  is  also  addressed  in 
Section  IV.  Conclusions  and  recommendations  are  presented  in  Section  V. 
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SECTION  II 

PLANE  PULSE  SCATTERING 

The  first  section  is  intended  to  be  a  brief  expos4  of  Lee's  theoret¬ 
ical  analysis  of  dilatation  wave  backscatter  from  an  object  insonified 
by  a  plane  dilatation  pulse  (Reference  7).  The  time-domain  first  Born 
approximation  is  presented  therein  to  linearize  the  resulting  expression 
for  the  scattered  displacement  field.  A  product  of  this  study  is  pre¬ 
sented  in  the  following  section  where  a  result  obtained  in  Reference  7 
for  scattering  from  a  void  is  generalized  for  any  homogeneous  anisotropic 
scatterer.  It  will  be  shown  that,  in  the  first  Born  approximation, 
certain  combinations  of  density  and  stiffness  of  the  host  and  scattering 
materials  render  the  scatterer  "transparent"  to  the  incoming  waves. 

1 .  THE  FIRST  BORN  APPROXIMATION 


Following  the  development  of  Lee  in  Reference  7,  consider  a  piece 
R  of  linear  elastic  material  bounded  by  the  surface  3R.  A  transducer 
located  on  a  subset  A  of  3R  launches  a  dilatation  pulse  which  encounters 
a  region  B  of  different  linear  elastic  material  located  somewhere  within 
R,  as  shown  in  Figure  1.  If  B  is  far  enough  away  from  A,  and  if  the 
radii  of  curvature  of  A  are  large  enough  compared  to  the  dominant  wave¬ 
lengths  of  the  pulse,  both  the  wave  incident  upon  B  and  the  scattered 
wave  received  by  the  transducer  can  be  approximated  by  plane  waves.  The 
analysis  which  follows  is  not  initially  restricted  to  backscatter,  i.e., 
the  observation  vector  x  may  point  in  any  direction.  (This  would,  of 
course,  require  another  tranducer  to  receive  the  scattered  waves.)  The 
condition  for  backscatter  is  imposed  later  in  the  section. 


The  material  densities  and  stiffness  tensors  are  given  by 
pO 


(x)  =|po 

1p°  + 


and 


cijkl(-}  = 


o 

°  +  A  p 

x  e  R-B 
*  — 

,  x  e  B 

(1) 

ijkl  5  *°«ij 

6kl  +  p°(6ik  6jl  +  6i1 

1  fijk^  -  e  R'B 

ijkl  +  ACijkl 

,  x  e  B 

!  _  % 

(2) 
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Figure  1.  Scattering  Configuration 
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where  X°  and  y°  are  the  Lame  constants  in  R-B  and  6..  is  the  Kronecker 

*  J 

delta  defined  by 


(3) 


The  incident  pulse  is  of  the  form 


UlnC(x,t)  =  V<j> 


e  .  x  +  B* 
t  -  - 

a  o 


(4) 


where  a0  is  the  dilatation  wave  speed  in  the  region  R-B  and  e  is  a  unit 
vector  characterizing  the  direction  of  propagation  of  the  pulse.  The 
constant  B  allows  the  time  origin  to  be  adjusted  so  that  t  =  0  corresponds 
to  the  instant  that  the  pulse  first  encounters  B.  Thus, 


B  =  -min(e  •  x)  (5) 

In  Equation  4,  the  scalar  potential  $(s)  satisfies 

<|>(s)  e  (-»,«>)  (6) 

and 

<f»(s)  =  o,  s  <  o  (7) 

A  minor  modification  of  Equation  5.4.2  of  Reference  7  defines  the 
pulse  length  of  Tp  as  that  time  after  which  a  specified  fraction  e  of 
the  total  energy  in  the  pulse  remains  according  to 

a>  2 

fj  Ms)  |  ds  =  e  £  1 
P 

The  displacement  field  then  satisfies 
inc  ,  . 

U  (X  ,t)  ,  -co  <  t  £  o 

inc  ,  sc  , 

u  (x,t)  +  u^  (x,t)  ,  t  >  o 


(8) 

(9) 


An  application  of  Love's  integral  identity  (Reference  8)  gives  the 
scattered  field  as 


(Ac..  usc  ),. 
klmn  m,n  1 


(10) 


where  (x,t)  is  a  body  force  field  which  is  dependent  upon  the  incident 
displacement  field  along  with  the  density  and  stiffness  perturbations 
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^(x.t) .  (bc..klU^),,  -  ipu;nc 


mi 


In  Equation  10,  the  dots  indicate  differentiation  with  respect  to  time, 
and  the  functional  U.^  is  given  by 


U,k(>.t|w(-)]  -=  3Zl\  _  «!k]  j/b”  aw(t .  ar)do 

L  r  rJ  l/a0 


(12) 


where 


r  =  |r]  =  |x  +  i\ 


(13) 


and  b0  is  the  shear  wave  speed  in  R-B.  Equation  10  shows  that  the 
scattered  field  arises  in  an  obviously  nonlinear  manner.  It  may  be 
possible  to  solve  for  ijsc(x,t)  by  an  iterative  solution  of  Equation  10; 
however,  a  simpler,  linear  problem  results  by  assuming  that  the  inter¬ 
action  of  the  scattered  field  and  its  derivatives  with  the  scatterer  is 
much  smaller  than  that  of  the  incident  field.  This  requires  that 


SC  ••  s  c 

Aci j^i »  Ap.  u.  (x,t),  and  u  (x,t)  be  small  in  some  sense  so  that 


Iff  U1k[t,t|  (ickli(nu“n)  -  lp'uskc]  dv) 

«MJ  1fk[r,t|bk(Z,.)]  dvy 


(14) 


Equation  14  is  a  statement  of  the  "time-domain  first  Born  approxima 
tion"  for  elastodynamic  scattering  which  then  gives  the  scattered  field 
as 


UtC  '  dvy 


(15) 
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Useful  results  may  be  obtained  by  considering  far-field  dilatation 
wave  backscatter  resulting  from  a  plane  dictation  pulse  input.  Generally, 
both  dilatation  and  shear  waves  will  be  scattered  due  to  mode  conversion 
at  the  scatterer's  boundaries;  however,  at  distances  which  are  large 
enough,  time-gating  may  be  employed  to  observe  only  the  response  due  to 
the  faster- travel ing  dilatation  waves. 

Be  expanding  Equation  15  with  Equation  12,  discarding  terms  which 

p 

are  0(x  )  and  terms  corresponding  to  shear  waves,  and  applying  Gauss' 

divergence  theorem  to  what  is  left,  the  far-field  scattered  dilatation 
field  uf(x,t)  is  obtained  as 


asc  -1 
u.  =  - 


*  vk  ApW  “k"^’1  -  i)] 


ixkxl  /u\*inc  /  .  £  \ 

a„  cklmn  ^um,n  [%-*  ~  a „J 


where 


x  =  x 


-  *  TIT 


Taylor's  expansion  of  r  and  order  of  magnitude  arguments  allow  further 
simplification  of  Equation  16.  The  resulting  expression  for  the 
scattered  dilatation  field  for  an  incident  plane  dilatation  pulse  is 


4irp0ao^C(x.t)  = 


*  A  A 

?  *„ 


x  e. 

(y)e  e  +  Ap(y) 
mn  ■*-  m  n  ae  ^ 


where 


y(t)  =  *(t) 


is  the  amplitude  of  the  incident  displacement  field,  and 


t'  =  t  - 
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2.  A  HOMOGENEOUS  ANISOTROPIC  SCATTERER 

The  result  obtained  in  Equation  19  is  a  statement  of  the  first  Born 
approximation  for  the  most  general  anisotropic  inhomogeneous  scatterer 
embedded  within  an  isotropic  homogeneous  host.  Reference  7  goes  on  to 
consider  further  details  for  the  special  case  of  scattering  from  a  void. 
In  this  report.  Equation  19  is  evaluated  for  a  homogeneous  anisotropic 
inclusion.  The  solution  has  the  same  form  as  that  for  the  void,  with 
the  addition  of  a  multiplicative  amplitude  factor  which  depends  upon  the 
materials'  elastic  constants. 


The  perturbations  attributable  to  the  scatterer  are  taken  as 


Ap  =  P  -  p' 


ACi jkl  =  Ci jkl  ‘  cijkl  (23) 

B  B 

where  p  and  the  c?^  are  constants  within  B,  and  c?^  is  as  defined 
in  Equation  2.  Substituting  these  into  Equation  19,  evaluating  with  the 
Kronecker  deltas,  and  setting  x  =  -e  for  backscatter,  yields 

4irp°aou!G  (x,t)  = 


T  W  p  -p  - 


B  „  X°  +  2y° 
-p  ’“2 


e,e,e  e  c. , 
k  1  m  n  klmn 


Consider  the  product  e in  Equation  24.  If  the  c^  ^  are 
contracted  according  to  the  method  described  by  Nye  (Reference  9),  the 
product  can  be  written  out  in  its  full  glory  (for  future  reference)  as 
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*ie j*kelCi jkl  “  e1c11  +  e2C22  +  *3C33 


*  4e24?e44  *  4eWc55  *  4eMC66 


4  *  2«?'3<=,3  4  2‘l'3e23 


4  4.2«jC15  4  4«2e2c,6  *  4.,.’c26 


4  4«’«jC24  4  *e2e3c34  *  4e1e3c35 


j  2  2 

+  4 e  e  e  c  +  4e  e  e  c  +  4e  e  e  c 

1  2  3  14  1  2  3  25  "l  2  3  36 


+  8e1 e2*lC4S  *  ee1*2e3C46  4  8e1*2*3C56 


(25) 


Now,  since 


.  >°  4 

°  '  a! 


and 


cf-|  “  P°a'o 


(26) 


(27) 


Equation  24  can  be  rewritten  as 


a;c(x.t)  = 


4ira 


dv 


(28) 


Equations  25  and  28  together  give  the  backscattered  dilatation  wave  ampli 
tude  from  any  homogeneous  anisotropic  scatterer  in  the  first  Born  approxi 
mation.  This  result  has  the  same  form  as  that  obtained  for  a  void 
(Reference  7),  i.e.,  the  backscattered  dilatation  wave  amplitude  is 
proportional  to 
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♦  SC(f)  = 


t'  - 


2 e  •  £  +  B 


a. 


dvy) 


(29) 


It  is  shown  in  Reference  7  that  performing  the  integration  results  in 

SC/*’'  =  Y  f  y(t'  -  t)  A(t)c1t  (30a) 


♦  (t’)  =  y  J  fU*  -  t)  A(t )dr 
0 


or 


t* 


4>SC( t 1 )  SY  f  Y^t*  "  A(f)di 


(30b' 


(2)  0  (2) 

for  both  y(t)  e  c  (-»,  «)  and  A(t)  e  c  (-»,  »),  where  Equation 


30a  is  obtained  for  y(0)  =  0  and  Equation  30b  for  A(0)  =  0.  To  obtain 
these  results,  a  coordinate  system  has  been  introduced  with  the  z  axis 
parallel  to  e  and  the  time  t  is  defined  by 


6  +  2z 


(31 


Note  that  t(z  .  )  =  0  so  that 
mi  n 


6  =  -2z, 


min 


(32 


The  function  A(t)  Is  a  function  which  describes  the  scatterer’s  cross- 
sectional  area  as  a  function  of  twice  the  one-way  transit  time  of  the 
plane  wave  passing  through  the  scatterer.  The  dependence  of  A(t)  upon 
twice  the  travel  time  makes  sense  physically,  since  the  backscattered 
portion  of  a  pulse  which  travels  some  distance  into  the  scatterer  must 
also  traverse  the  same  distance  back  through  the  scatterer. 


By  combining  Equations  28  and  29,  the  scattered  displacement  field 
can  be  expressed  as 


GiC(*.t)  - IT  M  ♦SC(t) 

4iraoX 


(33 


where  the  material -dependent  amplitude  M  is  defined  as 


B 


B  e.e,e_e_c.. 
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It  is  instructive  to  consider  a  homogeneous  isotropic  inclusion  B 
in  order  to  obtain  HfS0tr0piC  in  terms  of  familiar  engineering  quantities. 


It  is  convenient  to  define 


XB  =  X°  +  AX 
yB  =  y°  +  Ay 


E  =  E°  +  AE 


v  =  v°  +  Av 


I 


where  E  and  v  are  Young's  modulus  and  Poisson's  ratio,  respectively. 

g 

Substitution  in  Equation  34  of  an  equation  for  analagous  to  Equa¬ 

tion  2  for  isotropic  media  leads  to 

B  c?, 

"isotropic  -h  *  ^  '  2  (39> 

Note  that  this  expression  reduces  to  the  result  obtained  in  Reference  7 
for  scattering  from  a  void,  namely 


M  .  .  =  -2 
void 

Equation  39  is  alternately  expressed  in  terms  of  E  and  v  as 

M  .  oB  ,  EB(1  -  vR)(l  *  v°)(l  -  2v°)  2 

isotropic  p°  e°(1-  v°)(1  +  vB)(l  -  2vB) 


which  is  linear  in  AE  and  Ap 


Consider  the  ratio  in  Equation  41 


R  =  \  ■  -y  Ml  t  *  Ml  -  cv  ! 
v  "  (1  -v°)(l  +  vB)(l  -  2vB) 

g 

Substituting  for  v  from  Equation  38  and  expanding  yields 


n  -  v°-  av 
1  -  v° 


T^('  -W) 
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If  Av  satisfies 


Av  « 


1  +  v°|  ,  -1  <  v°  <  -1/4 


|  v°I  ,  -1/4  <  v°  <  1/2 


Rv  can  be  approximated  by  keeping  the  first  two  terms  in  the  geometric 
series  for  the  last  two  factors  in  Equation  43 

«  *  i  ♦  y°> —  4.  <45> 

(1  -  v"2)(l  -  2v“) 

Thus,  M  is  also  linear  in  Av  for  small  Av. 

The  possibility  of  a  scatterer  appearing  "transparent"  to  a  plane 
dilatation  pulse  in  the  first  Born  approximation  is  suggested  by  Equations 
34  and  39.  By  setting  M  =  0,  one  form  of  the  "transparency  condition" 
for  an  isotropic  scatterer  is  obtained  from  Equation  41  as 


’E®  (1  -  vB)(l  +  v°)(l  -  2v°)  p< 


M] 
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SECTION  III 

SCATTERING  FROM  SPHERES  AND  CYLINDERS 

The  results  of  the  previous  section  for  a  scatterer  of  arbitrary 
shape  are  applied  to  spheres  and  cylinders.  First,  the  cross-sectional 
areas  of  the  scatterers  as  a  function  of  time  are  determined.  A  model 
for  the  amplitude  of  an  Incident  displacement  pulse  as  a  function  of 
time  is  then  Introduced.  The  form  chosen  for  the  pulse  Is  not  only 
mathematically  attractive,  It  also  closely  models  the  pulse  which  Is 
produced  by  physical  transducers.  With  the  area  and  pulse  functions 
determined,  $sc(t)  Is  evaluated  according  to  Equation  30a  for  spherical 
and  cylindrical  inclusions  of  arbitrary  Isotropic  linear  elastic  constants. 
The  analytic  character  of  these  solutions  is  Investigated,  followed  by  a 
look  at  scattering  In  the  long  wavelength  limit. 

1 .  CROSS-SECTIONAL  SCATTERING  AREAS 


Consider  a  plane  dilatation  wave  pulse  traveling  In  the  positive  z 
direction  and  Incident  upon  a  spherical  scatterer  B,  of  radius  s,  as 

A 

shown  In  Figure  2a.  The  cross-sectional  area  A$(z)  "seen"  by  the  wave 
Is  the  area  of  the  family  of  circles  subtended  by  the  plane  on  the  sphere 
2  2 

ir(s  -  z  )  ,  -s  <  z  <  s 

0  ,  otherwise  (47) 


Recalling  Equations  5  and  31  for  B  and  t,  the  area  As(t)  Is  obtained  as 


,  0  < 


,  otherwise 


(48; 


A$(t)  Is  parabolic  as  Illustrated  in  Figure  2b. 

Now,  consider  a  dilatation  wave  pulse  incident  upon  a  cylinder,  where 
the  vector  e  Is  perpendicular  to  the  axis  of  symmetry  of  the  cylinder, 
as  shown  In  Figure  3a.  If  the  radius  of  the  cylinder  Is  c,  and  the  length 
of  the  cylinder  over  which  the  wave  can  be  approximated  by  a  plane  wave 
Is  L,  the  cross-sectional  area  Is  the  family  of  rectangles  described  by 


00 


0.  50 

a0T/4s 


1.  00 


b. 


Figure  2.  Plane  Wave  Interaction  with  a  Sphere 
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KM 


21  Jr2  -  z 


,  "c  i  2  I  c 
,  otherwise 


or 


AJt) 

c 


,  0  <  T  < 


4c 


,  otherwise 


(49) 


(50) 


which  describes  half  of  an  ellipse,  as  Illustrated  In  Figure  3b. 

SC 

Since  <|>  (t)  is  given  in  Equation  30a  as  a  convolution  of  a  pulse  am¬ 

plitude  function  with  A(t),  Figures  2b  and  3b  suggest  (by  considering  graph¬ 
ical  convolution)  that,  perhaps,  spheres  and  cylinders  will  give  similar 
scattered  time  waveforms.  This  topic  is  considered  in  more  depth 
In  the  section  on  long  wavelength  scattering.  However,  the  finite  slope 

of  the  sphere's  area  function  at  t  =  0  and  t  =  —  compared  to  the  infinite 

a  o 

slope  of  the  cylinder's  area  function  at  the  corresponding  points  will 
be  shown  to  give  characteristically  different  scattered  responses  when 
the  long  wavelength  limit  is  not  valid.  Note  that  it  is  precisely  the 
infinite  slope  Ac ( 0)  which  requires  the  use  of  Equation  30a  instead  of 
Equation  30b  to  obtain  4>sc(t)  for  the  cylinder. 


It  will  prove  convenient,  in  what  follows,  to  introduce  a  normalized 

time 


where 


t' 


_t 

2b 


2s 

ao 

,  sphere 

2c 

aQ 

,  cylinder 

(51) 


(52) 


This  allows  the  incident  pulse  to  be  characterized  in  terms  of  the  width 
of  the  scatterers'  cross-sectional  area  time  functions.  Note  that  2b 
Is  the  time  it  takes  the  wave  to  completely  traverse  the  scatterer  in 
the  incident  direction  and  again  in  the  backscattered  direction. 
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Using  Equations  51  and  52  In  Equations  48  and  50,  the  normalized 
area  functions  are 


As(t'  ) 


and 


Ac(t') 


2k2  l  . 

*a0t>  (t 


0 


2.  INCIDENT  PULSE  MODEL 


,  0  ^  T  '  <_  1 
,  otherwise 

,  0  <_  t  1  <_  1 
,  otherwise 


(53) 


(54) 


A  mathematically  tractable  and  physically  representative  model  of  the 
Incident  plane  dilatation  pulse  Is  given  as 


r(t)  = 


te"qtsinut 

0 


,  t  >  0 
,  t  <  0 


(55) 


where  w  Is  the  center  frequency  of  the  pulse.  The  factor  of  t  Is  used  to 
obtain  y(0)  =  0  for  use  of  Equation  30a.  Differentiating  twice  with 
respect  to  time 

(q^  -  u>^)te  ^slnwt  -  2qwte  ^coswt 


T<t)  * 


-2qe”^s1nwt  +  2we-c^coswt 


,  t  >  0 


'  0 


,  t  <  0 


(56) 


Using  Equations  51  and  52  In  Equations  55  and  56,  the  normalized  incident 
pulse  amplitude  functions  are 


r(f) 


2bt'e  sin  u'  t' 

0 


,  t*  >  0 
,  t'  <  0 


(57) 
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2b(q^  -  a>2)t'e  sinui't'-Abqwt'e"2^  cosw't' 
-2qe"^^c^t  sinaj't'  +  2we  cosw't',  t'  >  0 

0  ,  t'  <  0 

(58) 

where 


a) 1  *  2bw 


is  the  normalized  radian  frequency.  In  terms  of  the  period  T  =  -  ir 

flV  * 


T.  -  L_  (60) 

1  "  2b 

gives  the  period  of  the  pulse  relative  to  the  time  width  of  the  scatter- 
ers*  area  functions. 


It  is  revealing  to  put  this  time  domain  pulse  normalization  in  the 
perspective  of  the  frequency  domain  wavenumber  normalization  ka  employed 
throughout  the  literature  where  a  is  the  characteristic  dimension  of 
the  scatterer,  and  k  =  ~  is  the  wavenumber.  If  a  corresponds  to  either 
radius,  s  or  c,  recalling  Equations  52  and  6(T  for  b  and  T'  yields 


W,  a  -JL_ 

m  ZV  (61) 

3.  SOLUTION  FOR  THE  SPHERE 

With  all  the  pieces  now  available,  the  normalized  solution  for  the 
backscattered  amplitude  function  4>s  (t*)  for  scattering  from  a  sphere 
may  be  obtained.  Either  form  of  Equation  30  may  be  used  to  obtain 
4>s  (t')  since  As(0)  =  0.  For  the  form  given  by  Equation  30b,  the 
required  second  derivative  of  the  area  function  is 

ira?b2[6(T')  -  2  +  6(t'  -  D]  ,  0  <  t'<  1 

A  (f)  = 

0  ,  otherwise  (62) 
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The  delta  functions  greatly  simplify  the  convolution  with  the  incident 

field  in  Equation  30b  and  only  the  convolution  with  the  constant  (-2) 

must  be  worked  out  in  detail.  The  situation  is  not  so  pleasant  for  the 

cylinder  and  Equation  30a  is  required  due  to  the  singular  behavior  of 
• 

Ac(t')  at  t'  =  0  and  at  t*  »  1.  In  order  to  provide  a  unified  treatment 
of  scattering  from  both  objects,  the  form  of  Equation  30a  is  chosen, 
obtaining  the  solution  in  terms  of  the  area  function  convolved  with 
second  derivative  of  the  incident  pulse. 


Substituting  Equations  53  and  58  into  Equation  30a  produces 

^C(t')  =  2traob3|‘|>(q2  -  o>2)(f  -  T*)e"2bq(t'  '  T,)sW(t'  -  x‘) 

V  -2bqai(t'  -  r')e'2bq(t>  '  ^ costa 1  (t 1  -  x') 

-qe"2bq(t’  ‘  T*  ^sinu)'  (t1  -  x')  +  ^‘^^‘"'^cosa)'  (t1  -  x')}; 
x  tx'  -  x'2  }  dx'  ,  x'  >  0  (63) 


where  dx  *  2bdx'  is  employed,  and 
(f  ,  0  <  f  <  1 

(l  ,  f  >  1  (64) 

which  limits  the  Integration  to  that  portion  of  the  sphere  whichis  en¬ 
countered  by  the  wave  at  time  t'.  Introducing  the  complex  exponent 

K*  =  2bq  ±  iw'  (65) 

and  using  Euler's  formulae  for  the  sine  and  cosine  functions.  Equation 
63  is  expanded  and  rearranged  to  obtain 


♦ie(t;) 

2na^b2 


-i{A*(t')  -  B*(t')}.?  T'eK"Tldx’ 

+i|A(t')  -  B(t* )J  -T  x ' e*+T ' dx ' 

+^{A*(t')[f  +  1]  -  B*(t')}^  x'2  eK_T  dx' 
-i|A(t')[t'  +  1]  -  B(t')|  }  x'2eK+Tdx' 
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-i  A*(t' )  ^*r'“eK  T  dx' 
+i  A( t *  )^r’3eK  T  df 


,  t'  >  0 


where 


«»•>■  +  fb2qu]  e~K+t' 


(66) 


(67) 


and 


B(t') 


■[ 


m  *4]  e-"+t' 


(68) 


The  asterisk  denotes  complex  conjugation  and  l  =  *^T.  The  integrals  in 
Equation  66  are  easily  handled  by  defining 


C .(a)  =  /  t‘V'  1  dx* 


°  J.RV 

'j'  *  o  T 

Integration  of  Equation  69  by  parts  leads  to  the  recursion  relation 

aJ@K+a  -  cj  i(a) 

K+  J-1 


(69) 


c  (a)  = 

J  k+ 


(70) 


which  is  particularly  useful,  since  C-j(a)  is  readily  evaluated.  Note 
also  that 


C*(a)  =  ./x'^e^  T  dt' 
J  o 


(71) 


so  that  Equation  66  can  be  rewritten  as 
~sc(t' ) 

--y?  3  4A(t')C,(a)  -  A*(t')C*(a)l  t* 

2Tra„b^  L  1  1  J 


'  / 


[b ( t '  )C]  (a)  -  B*(t')C*(a)] 
-t[A(f)C2(a)  -  A*(t')C*(a)][r  +  l] 
+i[B(t')C2(a)  -  B*(t')C*(a)] 

)C3(a)  -  A*(t')C*(a)] 


,  t'  >  0 


(72) 


20 


AFWAL-TR-82-4044 


If  Im{*)  denotes  the  Imaginary  part  of  a  quantity,  the  backscattered 
wave  amplitude  from  a  sphere  is  finally  obtained  as 


2  |  Im[A(t' )C] (a)]  f 
- Im[B ( t 1 )C] (a)] 

-Im[A(f)C2(a)]  [f  +  1] 

+Im[B(t')C2(a)3 

-Im[A(t')C3(a)]  J  ,  f  >  0  (73) 

where,  recall,  a  equals  t'  or  1  according  to  Equation  64. 

While  Equation  73  is  not  particularly  revealing,  it  is  a  beautifully 
simple  result  that  is,  together  with  Equations  69  and  70,  quite  amenable 
to  numerical  evaluation.  It  is  easy  enough  to  observe  from  Equations  67 
through  70  and  Equation  73  that  <f>  ( t ' )  has  the  same  center  frequency 

as  y(t')  and  falls  off  as  t'e"^b^'  for  large  t'. 

4.  SOLUTION  FOR  THE  CYLINDER 


By  a  similar  development  as  that  for  the  sphere,  ( t ’ )  is  obtained 
for  the  cylinder  by  substituting  Equations  54  and  58  into  Equation  30a 
producing,  upon  rearrangement 


4Laob 


|[b2(q2  -  w2)Vy  -  bqjsinu't'  - |j2b2qt ' -bj  wcosu't’l 

x/e+2bqT'  /T*-T cosw’x’dx' 
o 

-|Jb2(q2-a)2)t'-b^cosa)'t'  +^b2qt 1  -bju>sina> ' t '  ^  e"2bqt 

x/  e+2bqT  A '-t ' 2  sinw'x'dx' 
o 

-  |b2(q2-u)2)sinu>'t'  -  2b2q  cosw't'  |  e"2^qt 
x/  e+2bqxf  T »/T ' _T •  2  cosw'x'dx' 


-2bqt' 
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+  j[ b2(q2-<o2)cosa>* t*  +2b2Mqs1n«'t'  j  e'2bqt' 

xJr-  e+2bqx ' T i ^ i _T sinu'x'dx1  ,  t'  >  0  (74) 


t/ 


where,  again  a  =  t'  or  1  according  to  Equation  64.  The  integrals  in 

*■  jr 

Equation  74  lead  to  a  solution  for  $  (t')  which  is  not  as  straightfor- 

~  sc  ^ 

ward  as  <j»s  (t')  obtained  for  the  sphere.  Consider  the  Integral 

,  r 

(f )  =  ^  e2bqT'*/rr  (I-t1)*5  cosoj't'  dr'  (75) 

Expansion  of  the  radical  in  a  binomial  series  gives 

M*’>  =  £  [e2bqVIVcosu’.t'  £  (-DJ(  j)*'1]*1  (76) 

provided  0  <_  t'  <•  1.  Now,  since  the  series  Is  uniformly  convergent  with¬ 
in  the  interval,  the  summation  and  Integration  may  be  interchanged  to 
obtain 


) 


(77) 


where 


\c.  (t')  h  ^  T,is+^  e2bqT  cosu't'  dx' 
J  ^ 


(78) 


Similarly, 


I2(f)  =  ^  e2btlT '/^r*  (i.T-)14  sinw'x'  dx' 


(79) 


1/ 


becomes 


I2(f) 


,4  '-"‘(s)*; 


(t1) 


where 


(80) 


A3(f) 


•  X,>s+j  e2bqTJ  sinw’x’  dx’ 


(81) 
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The  functions  A ^ ( t * )  and  ( t ' )  lend  themselves  to  recursion  relations 
which  make  them  computationally  attractive.  Integrating  Equations  78  and 
81  by  parts  yields 


c  t|Js+j  e2bqt’  r  l 

A^(t')  = - r — 5 —  [2bqcosu)'  t'  +  aj'slnoj't'] 

J  (2bq)z+u,'z 


0s+.i)2bq 
(2bq)Z+a>' 2 


■  15  V  2Aj-l(t>) 

(2bqr+w,Z  J"1 


(82a) 


and 


e2bqf 

(2bq)Z+  a.'2 
.  (%*J)2bq 
(2bq)2-Ku*  2 


[2bqs1nu't' 


w’coSw't*] 


k  *  j 

- 9 — 7 

(2bq)Z+ur 


*j.it**> 


(82b) 


The  appearance  of  the  recursion  Index  j  In  the  numerators  of  the  forward 
recursion  re1at1ons((Equat1on  82)  Is  not  desirable  since  truncation  errors 
Introduced  In  the  numerical  evaluation  of  the  functions  will  be  amplified 
as  j  gets  large.  A  numerically  stable  evaluation  scheme  results  from 
manipulation  of  Equation  82  which  leads  to  the  reverse  recursion  relations 


Aj(t,)  =  1?l+j[t,3/2+J  e2bqt'  C0SW,t'  +u'Aj*l(t,>  “2bqAj+1  (t‘ )]  (83a) 
Aj(t,)  "  37kj[t,3/2+j  e2bqt'  s1nw,t'  +u'A5+i(t,)  ~2bqAj+1  ( t ' )]  (83b) 
Note  that 

11m  A^(t')  *  11m  A*(t')  -  0  (84) 
j  -*■  »  J 

which  Is  easily  seen  by  considering  Equations  78  and  81  remembering  that 
0  <  t'  <  1. 

Because  of  Equation  84,  the  reverse  recursion  relations  are  also 
convenient  In  that  an  exact  expression  Is  not  needed  for  one  of  the 
recurring  functions  In  order  to  start  the  recursion.  For  large  j. 


3/2+j+2bq 


j  - 


(86b) 


It  would  be  nice  to  somehow  Insure  that  the  functions  are  being 
generated  correctly  as  j  Is  decreased  from  Its  starting  value  to  zero. 

A  great  deal  of  confidence  could  be  placed  In  the  scheme  If  the  final 
recursions  to  Ag(t')  and  aqU')  matched  a  known  "correct"  answer  to 
within  some  desired  degree  of  accuracy.  Surprisingly  enough,  the  zeroeth 
Integrals  can  be  evaluated  In  closed  form.  This  requires  a  momentary 
diversion  to  consider  the  Integral 

I  =  £eks  s**  ds  (87) 

2 

Making  the  substitution  s  =  u 

I  ■  2  u2  eku2  du  (88) 


which  can  be  written  as 


Furthermore,  the  substitution  v  *  in,  where  i  *  renders 
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which  Is  recognized  to  contain  a  form  (Reference  10)  of  the  complex  error 
function 


(92) 

Performing  the  Indicated  differentiation  (Reference  10) 


1  =  T  ekt  +  erf  (93) 

Back  to  evaluating  A*(t')  and  A{j(f),  Euler's  formulae  for  the  sine 
and  cosine  functions  In  Equations  78  and  81,  and  an  application  of  Equa¬ 
tion  93  to  both  yield 


c  /£'  p29bt’ 

AQ(t' )  =  ^ - v  [2bqcoso)'t'  +  wslnu't'j 

(2bq)  +  to ' 

^  f erf  (i/k+t1)  erf 

4  I  (K+)3/2  '  (K-)3/2  J 

and 


(94a) 


/F  e2bqt> 

(2bq)2+to'2 


[2bqs1nw't' 


ojcosoj  '  t  'J 


+ 


A 

4 


erf  (t/K^t* ) 
(K+)3/2 


erf  (i/Ft1 ) 

~7n3/2  . 


(94b) 


where,  as  before  K±  =  2bq  ±  iu\  Equations  74,  78,  and  81  can  now  be 
combined  to  give  the  solution  valid  for  0  <  t'  <  1,  l.e.,  when  the  front 
edge  of  the  pulse  Is  within  the  scatterer.  The  result  Is 


4taob 


=  e 


2bqt' 


r  (-dj 

j*0 


(Dr*, 


(t')  A^(t') 

J 


-  A2(f)  A*(f)-  *3(f)  Aj+1  (t' )  +  A4(f)  A*+1(t*)]_ 

0  <  t'  <  1  (95) 


where 
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A]  ( t * )  »[b2(q2-w2)  t'  -bqjslnw't*  -  [2b2qu,t'  -  bwj 

J  cosui  *  t  *  ( 96a ) 

A2(t')  *[b2(q2-u>2)  t'  -  bqjcosw't'  +^2b2qu>t‘  -  buj 

slnw't'  (96b) 

A3(t')  =  b2(q2-&>2)  sinui't'  -  2b2qu>tosa>'t' 

(96c) 

A^(t’)  =  b2(q2-a>2)  cosu't'  +  2b^qtu6lnoj't' 

(96d) 

The  functions  A ^ ( t * )  and  ( t * )  are  started  by  Equation  86,  evaluated 
by  Equation  83,  and  compared  to  Equation  94  to  Insure  a  desired  degree 

of  accuracy. 

The  solution  for  t'  ^  1  must  still  be  obtained.  Equation  74  for 

SC 

i  ( t ' )  Is  equivalent  to  a  form  more  Indicative  of  the  convolution  opera- 

C' 

tlon  (Indicated  by  an  asterisk),  and  thus,  amenable  to  Laplace  transform 


methods . 

The  desired  form  of^c(t')  Is 

SC(t,) 

4La2b 

■  b2(q2-a>2)f1(t')  -  2b2qu>f2(t')  -  bqf3(f)  +  ba>f4(f) 

(97) 

where 

f , ( t ' )  =  (t'  e‘2b<,t,s1nu,t')  *  e(t') 

(98a) 

f2(t')  -  (t'  e’^'cosw't')  *  e(t') 

(98b) 

f3(t')  =  ( e~ 2 b<qt ' s  1  n<*> ' t ' )  *  0(t') 

(98c) 

f4(f)  =  (e‘2bqt,cosu,t')  *  Q(t’) 

(98d) 

and  where 

/  /t'  -  t'2  ,  0  <  t'  <  1 

-e(t')  ■  j 


0 


,  otherwise 


(99) 
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If  L  (•)  denotes  the  or.e-slded  Laplace  transform  of  a  quantity,  the 
required  transforms  are  evaluated  In  terms  of  the  transform  variable  s  as 


L  [t  e"qt  sin  wtj 


L  [e_<^  cosutj 

and  from  Reference  11, 
L 


2b)(s  ♦  q) 
(s+q)2  +  w  1 


L  |"t  e~qt  cosutl 

_ 

(s+q)2  -  o>2 

L  J 

[(s+q)2  +  0,2] 

L  fe"qt  slnutl 

s 

0, 

(s+q)^  +  u/ 
s  +  q 

(s+q)Z  +  u,2 


v/2  e“s/2  I.(s/2) 


(100a) 

(100b) 

(100c) 

(lOOd) 

(101) 


where  Iv(*)  Is  the  modified  Bessel  function  of  the  first  kind  of  order 
v.  Denoting  L[fj(t')]  by  F^(s),  for  i  =  1,  2,  3,  4,  the  F^(s)  are  ob¬ 
tained  by  the  multiplication  of  the  appropriate  transforms  In  Equation 
100  by  the  transform  In  Equation  101.  The  f . (t ’ )  are  conveniently 
obtained  by  summing  the  residues  at  the  poles  of  e11,  F.(s),  which  poles 
are.  In  all  cases 


s1  =  -K1  =  -[2bq  ±  u] 

II  (fO 

Note  that  -  Is  an  entire  function  of  s  and,  therfore,  s  = 

s 

pole  of  the  F^(s).  Saving  the  details,  the  result,  valid  for 
Is  obtained  as: 


(102) 

0  is  not  a 
t'  >  1, 


^  *  f{[bq-b2  (q2-o,2)f]  Im  [  E(f  )] 

4LaZb  2  lL  J 

+  bo,  (1  -  2bqt' )  Re  [E(t’)] 
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b2  (q2-w2)  Im  [0(fj| 
2b2qw  Re  [®(t')]J 


,  t*  >  1 


where 


•  >  *  '■Kt(t'J,){  ^  [Vt>  *  <T)]( 


0  03) 


(104) 


Elf)  5  e-K+(t' -*5)  V|i 


(105) 


and  where  Re(*)  denotes  the  real  part  of  of  the  quantity. 

Equations  95  and  103  together  give  the  solution  for  scattering  from 
a  cylinder  for  all  t' ,  Including  the  critical  point  where  the  solutions 
are  pieced  together  at  t'  *1.  Again,  the  form  of  the  solutions  Is  not 
particularly  revealing,  however,  the  same  observations  as  those  made  for 
the  scattering  from  the  sphere  can  be  made  with  respect  to  the  scattered 
center  frequency  and  the  form  of  the  solution  for  large  t'. 

Given  the  simplicity  of  the  solution  for  t *  1  compared  to  the 

solution  for  0  <_  t'  <1,  one  Is  tempted  to  look  for  a  Laplace  transform 
solution  valid  for  all  t'.  However,  the  above  inversions  are  valid  only 
for  t'  >.1.  Consider,  for  example,  est  F,(s)  in  the  form 
‘  .  .  ....  I,(s/2>3 


St'  _  /  X  WIT 

e  Fj(s)  -  -T 


sls-sjxs  -  .-) 


(106)' 


The  asymptotic  expansion  of  Iv(z)  (Reference  10)  for  large  z  and  |arg(z)|< 
along  with  the  fact  that  I-j(z)  Is  an  odd  function  of  z  yields 


Il(s/2)|s| 


,  |  <  arg(s)  <  3  j 


(107) 


28 


IM|=I 


AFWAL-TR-82-4044 


Thus, 


WIT 

z7m 


e(f-l)s 

s(s-sJMs-s-) 


l<arg(s)<3|  (108) 


St1 

Since  the  poles  of  e  F_ ( s )  are  In  the  left  half  of  the  complex  plane, 
e  Fg(s)  required  to  vanish  as  |s|  ->  •  In  order  to  sum  residues  to 
obtain  f.(t').  However,  for  s  In  the  left  half  plane  (with  a  negative 
real  part),  i„i  11m  e  F,(s)  ■'«,  unless  t'  >  1,  which  therfore  limits 

|  S|  -*  »  J  — 

the  Laplace  transform  solution  to  the  region  t'  >_  1. 


5.  SMOOTHNESS  OF  THE  SOLUTIONS 

The  Interaction  of  the  Incident  pulse  with  scatterers  having  abrupt 
onset  and  terminating  boundaries  raises  questions  about  the  smoothness 
of  the  scattered  waveforms.  Recall  the  scattered  field's  second  time 
derivative  must  be  sufficiently  Integrable  so  that  Equation  14  Is 
satisfied  for  the  first  Born  approximation  to  be  valid.  The  solution 
for  the  cylinder  Is  particularly  suspect  due  to  the  Infinite  slope  of  the 
area  function  A  (t')  presented  to  the  wave  at  t’  =0  and  T'  =  1.  Thus, 

"SC  ^ 

<t>c  (t* )  Is  treated  here  first.  In  the  form  presented  by  Equations  97  and 
98. 

By  writing  out  the  convolutions  as  Integrals,  keeping  the  first  terms 
In  the  expansions  for  all  factors  In  the  Integrands,  and  then  performing 
the  Integrations,  the  result  obtained  Is 


Thus, 


and 


*c  (t  >  :  2bu>  t,3/2 

— 2 -  1“  1 

4Latb  6 


,  t'  «  1 


(109) 


*cC(t’>-  % 

-£-? -  -  bwt'^ 

4Latb 


,  t'  «  1 


(110) 


4Latb  T  1 


,  t'  «  1 


(111) 
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Note  that  4>sc(t'),  { t ' ) ,  and  <^C(t')  are  all  zero  If  t'  <  0.  Thus, 

«cr  C  C  C 

only  ♦  (t‘)  is  discontinuous  at  the  onset  of  scattering,  with  an  infi- 

c 

nite  discontinuity  from  a  t  singularity. 

Recall,  in  the  convolutions  of  Equation  97,  the  upper  limit  of  inte¬ 
gration  is  t'  for  0  £  t'  £  1.  Differentiating  with  respect  to  t'  accord¬ 
ing  to  Leibnitz'  rule  (Reference  12)  yields 


v.  ;  19  9  2  1  .  . 

-  =  -2bq  b^(q  -w  )  -  bwu>'  f,(t') 

4La„b  L  1 

+  [u>'b2(q2-u>2)  +  4b3q2a?]f2(t') 

+  [b2(q2-aj2)  +  2b2q2  -  bww']f3(t') 

-  £  4b2qu>  +  bqw  Jf.(t') 

+  bw/t'  -  t2  ,  0  £  t'  £  1  (112) 


where  the  f . ( t ' )  obviously  give  back  combinations  of  the  f.(t')  and  the 
last  term  results  from  differentiating  f^(t')  ala  Leibnitz.  For  t'  £  1, 


♦"(f> 


=  -2bq[b2(q2-u2)  -  b™']  g^t')  +[«i'b2(q2-«>2)  +  4b3q2<o]g2(t' ) 


[ 


9  9  9  2  9 

+  1  b  (q  -w  )  +  2b  q  -  bwui1 


]g3(t‘  )  -  [  4b2qw  +  bqw'] 


g4(t') 
(113) 


where  the  g^(t')  are  the  same  convolution  integrals  as  the  f. ( t ' ) ,  except 
that  the  upper  limit  of  integration  is  1,  instead  of  t'.  Since  g^(l)  = 
f . ( 1 ) ,  for  i  =  1,  2,  3,  4,  comparison  of  Equations  112  and  113  implies 

that  <L  ( t ’ )  is  continuous  at  t'  =1.  Moreover,  since  the  f . ( t ' )  e  Cv  ;(-~ 

c  ~sc  ill 

the  important,  physically  expected,  result  that,  $  ( t ' )  e  C'  '(-<*>,»)  is 

-er  ^ 

obtained.  Now,  differentiating  <j>c  ( t ' )  introduces  a  factor  of 

(t'  -t'2)-1^2  for  0  <  t'  <  1  which  obviously  means  that  <jisc(t')  is  not 

—  c 

continuous  at  t'  =  1.  In  fact,  .. 
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Similarly,  for  t'  >  1 


*!C(f) 

11m 

t,  +  1+  4bLaQ 


(115) 


Thus,  $*c(t')  has  an  infinite  discontinuity  (at  t' 
t'  =  0. 


1  just  as  it  does  at 


It  is  reasonable  to  question  the  validity  of  the  first  Born  approxi- 

tj  er 

mation  in  the  light  of  these  results  for  <t>c  (t')»  since  this  relates 
directly  back,  through  Equation  33,  to  Equation  14.  The  approximation 
is  not  threatened,  however,  since  the  t'-1*  singularity  is  sufficiently 
integrable  in  Equation  12  and  volume- integrable  to  insure  that  Equation 
14  holds. 


The  same  analysis  applied  to  the  sphere  yields  that 


♦|C(t*  )255  ^  t'2 


,  t'  «  1 


J|c(t')2;  but' 


,  t'  «  1 


(116) 


$*c(t') »  bu> 


,  t'  «  1 


Again,  the  second  derivative  is  discontinuous  at  the  onset  of  scattering 
although  finite  for  the  sphere.  The  discontinuity  bu>  is  alternately 
expressed  as 


bw  =  2  ka 


(117; 


where  ka  is  as  mentioned  in  the  section  on  the  incident  pulse  character¬ 
istics.  The  second  derivative  is  found  to  be  discontinuous  at  t'  *  1 
by  -bw,  and  more  importantly,  $sc(t’)  e  (-”,”)• 
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6.  LONG  WAVELENGTH  LIMIT 

The  time-domain  first  Born  approximation  can  already  be  thought  of  as 
a  long  wavelength  approximation  in  its  own  right.  This  seems  reasonable 
since  it's  an  approximation  based  upon  weak  interaction  of  the  incident 
field  with  the  scatterer;  weak  interaction  is  intuitively  appealing  for 
wavelengths  long  compared  to  the  size  of  the  scatterer.  However,  it 
appears  that  the  approximation  as  stated  in  Equation  14  could  be  satis¬ 
fied  by  sufficiently  small  perturbations  Ap  and  ac^i  even  in  shorter 
(or  at  least  resonant)  wavelength  regimes.  This  certainly  is  an  item 
begging  for  further  attention. 

The  long  wavelength  limit  is  stated  according  to  T*  »  1 ,  or  equiva¬ 
lently,  ka  «  In  this  case.  Equation  30a  is  particularly  revealing 
since  both  As(t)  and  A^t)  now  sample  the  second  derivative  of  the 
incident  pulse  over  sufficiently  small  intervals  to  produce  results 
similar  to  convolutions  with  a  delta  function.  Thus,  in  the  long  wave¬ 
length  regime,  spheres  and  cylinders  produce  essentially  identical 
backscattered  fields  which  have  pulse  amplitudes  given,  very  nearly,  by 
the  second  derivative  of  the  input  pulse  amplitude.  This  result  is 
obtained  independent  of  the  shape  of  any  small  scattering  object.  The 
most  significant  difference  between  the  backscattered  fields,  in  the  long 
wavelength  limit,  will  be  the  amplitudes  of  the  scattered  fields  due  to 
the  "strength"  of  the  delta  function.  This  strength  depends  upon  the 

size  of  the  object  and  is  manifested,  in  the  case  of  cylinders  and  spheres, 

2  3 

in  the  volumes  Lttc  and  4ns  /3.  Of  course,  the  material  dependent  ampli¬ 
tude  factor  M  introduced  earlier  plays  a  rrle  in  the  scattered  amplitude. 

As  an  example  of  how  the  delta  function  strength  anu  affect  a 
backscattered  response,  consider  the  problem  of  distinguishir  the  responses 
from  a  long  cylindrical  Inclusion  and  a  spherical  void  In  the  long  wave¬ 
length  limit.  The  scattered  time  functions  will  have  indistinguishably 
similar  form,  to  within  a  minus  sign  attributable  to  M,  and  only  the 
relative  amplitudes  may  be  strikingly  different.  The  cylindrical  inclu¬ 
sion  may  very  well  have  an  M  which  is  quite  small  compared  to  that  for 
the  spherical  void.  However,  if  the  length  L  over  which  the  cylinder  is 
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illuminated  by  the  approximately  planar  wave  is  large  enough,  the  response 
from  the  weakly  scattering  inclusion  can  dwarf  the  response  from  the 
relatively  strongly  scattering  void,  due  to  the  strength  Uc  which  is 
unbounded  in  L.  This  problem  is  of  considerable  interest  to  the  Air 
Force  for  non-destructive  inspection  of  fiber-reinforced  composite 
materials,  where  detection  of  gas  bubbles  (porosity)  in  the  epoxy  matrix 
(host)  is  desirable.  It  would  be  nice  if  the  cylindrical  fibers  were 
relatively  transparent  to  the  incident  dilatation  pulse,  as  compared  to 
a  spherical  void,  allowing  the  response  from  the  porosity  to  stand  out 
against  a  "background"  of  cylinders'  response.  Further  attention  is 
given  to  this  item  in  the  following  section. 
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SECTION  IV 
NUMERICAL  RESULTS 

Some  of  the  numerical  techniques  used  for  computer  evaluation  of  the 
scattered  responses  from  cylindrical  and  spherical  inclusions  are 
presented.  Responses  obtained  for  incident  pulses  having  various  normal¬ 
ized  center  frequencies  ka  are  then  presented  to  illustrate  the  differ¬ 
ences  and  similarities  between  scattering  from  cylinders  and  spheres. 

An  experimental  result  for  scattering  from  a  cylinder  is  presented  and  a 
comparison  made  to  the  response  predicted  by  the  theory.  Finally,  the 
transparency  of  the  fibers  in  a  fiber-reinforced  composite  material  is 
addressed. 

1.  NUMERICAL  TECHNIQUES 

Recall,  the  incident  pulse  amplitude  is  modeled  as  an  exponentially 
damped  sinusoid  according  to  Equation  55.  For  the  following  results, 
the  length  of  the  pulse  Tp  defined  by  Equation  8  was  chosen  to  be  the 
time  after  which  ten  percent  of  the  energy  in  the  pulse  remains.  With 
a  desired  pulse  length  and  center  frequency  specified,  substitution  of 
Equation  55  into  Equation  8  leads  to  a  transcendental  equation  for  the 
damping  q  which  was  solved,  most  conveniently,  by  the  simple  method  of 
bisection,  as  described  in  Reference  13. 

The  error  function  and  modified  Bessel  functions  of  complex  arguments 
were  evaluated  by  routines  based  upon  the  forms  of  the  functions  given 
in  Reference  10.  Power  series  and  continued  fractions  were  utilized  to 
provide  at  least  eight  significant  figures.  Asymptotic  expansions  were 
employed  in  the  appropriate  regions.  The  reverse  recursion  relations 
(Equation  83)  were  considered  to  provide  accurate  results  if  they  matched 
the  closed  form  expressions  (Equation  94)  to  seven  figures.  Finally, 

251  points  were  calculated  for  each  theoretical  plot  displayed  from  here  on. 

2.  TIME-DOMAIN  BACKSCATTER  WAVEFORMS 

A  typical  incident  pulse  and  its  second  derivative  are  shown  in  Figure 
4,  illustrating  the  particular  case  when  the  pulse  length  is  equal  to  the 


Figure  4.  Incident  Pulse  Waveforms 
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period.  This  may,  of  course,  be  specified  for  a  pulse  of  any  period. 

The  specific  case  of  scattering  when  T1  =  0.5  ( ka  -  tt )  Is  shown  in 
Figure  5.  The  most  interesting  difference  between  the  responses  from 
the  cylinder  and  sphere  occurs  around  the  point  t'  =  1.  It  is  at  this 
time  that  the  portion  of  the  wave  which  makes  the  round  trip  through  the 
entire  scatterer  is  felt  at  the  observation  point  where  the  wave  first 
enters  the  scatterer.  The  more  interestingly  "structured"  response  from 
the  cylinder  Is  due  to  the  more  abrupt  nature  of  the  change  in  density 
and  stiffness  which  the  wave  sees  as  a  result  of  the  infinite  slopes  of 
the  cylinder's  area  function.  While  the  first  Born  approximation  is  not 
guaranteed  (or  expected)  to  be  valid  when  ka  is  this  large,  these  results 
are  interesting  in  that  they  predict  marked  differences  in  the  backscatter 
from  cylinders  and  spheres  at  shorter  wavelengths.  Figure  6  illustrates 
the  scattering  when  ka  =  ^  ,  for  the  same  pulse  shape  (T  =  T ' )  as  that 
In  Figure  4.  These  results  illustrate  less  structure  around  t'  =  1  than 
for  ka  -  tt.  This  is  to  be  expected  since  increasing  Tp  and  T'  have  the 
effect  of  narrowing  the  bandwidth  of  the  frequency  spectrum  of  the  inci¬ 
dent  pulse  and  centering  it  around  a  lower  center  frequency.  This  effect 
then  carries  through  to  the  spectrum  of  the  scattered  pulse  and  is  evi¬ 
denced  in  a  time-domain  pulse  of  less  structure. 

3.  COMPARISON  TO  EXPERIMENTAL  RESULT 

While  the  first  Born  approximation  may  be  suspect  for  values  of  ka 
as  large  as  in  the  examples  above,  excellent  agreement  between  an  experi¬ 
mental  result  and  theoretical  predictions  for  smaller  ka  will  be  presented. 

The  experiment  conducted  at  the  Materials  Laboratory  is  illustrated 
In  Figure  7.  A  voltage  waveform  f(t),  input  to  the  transducer,  pro¬ 
duces  a  plane  dilatation  pulse  which  Interacts  with  a  cylindrical  void 
(radius  50.8  ym)  located  inside  a  piece  of  linear,  homogeneous ,ancj  iso¬ 
tropic  elastic  material.  The  scattered  wave  gives  rise  to  an  output 

voltage  g  (t).  The  Incident  wave  also  reflects  from  the  planar  "back 
sc 

wall"  producing  the  response  g  (t).  This  response  can  be  considered 

Oh 

to  be  unperturbed  by  the  Incident  and  reflected  wave  traveling  through 
the  scatterer,  since  the  magnitude  of  the  scattered  field  is  much  smaller 
than  the  Incident  and  reflected  fields.  It  is  shown  in  Reference  14 
that  treating 
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Figure  5.  Backscatter  for  ka 


Figure  6 


b. 

Backscatter  for  ka  *  it/ 2 
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Scattering  Experiment 
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9sc(t)  =  9bw{t)  *  A(t)  (118) 

is  completely  equivalent  to  treating  the  convolution  (Equation  30a)  for 
the  scattered  field. 


Figure  8  shows  agreement  between  theory  and  experiment  for  scattering 
from  the  cylindrical  void  when  ka  =  0.32  (T'  =  4.85).  The  solid  curve 
in  Figure  8a  is  the  back  wall  reflection,  modeled  numerically  by  the 
broken  curve  given  by 

y (t)  =  t(e"qt  -e'pt)sino)t  (119) 

The  length  of  the  abscissa  is  1  ysec  for  all  graphs. 

The  scattered  response  and  the  numerical  prediction  are  similarly  illus¬ 
trated  in  Figure  8b.  Such  agreement  between  theory  and  experiment  lends 
credence  to  this  approach  to  elastic  wave  scattering.  The  lower  center 
frequency  of  the  reflected  pulse  can  be  attributed  to  high  frequency 
attenuation  present  in  the  host  material  and  not  accounted  for  in  the 
model.  The  scattered  pulse  is  not  affected  as  much,  since  it  travels 
through  only  about  one-half  the  distance  that  the  reflected  pulse  travels 
through.  Figure  8c  shows  the  second  time  derivative  of  Equation  119. 
Figures  b  and  c  Illustrate  the  result  predicted  earlier  that  the  scattered 
field,  in  the  long  wavelength  regime,  is  very  nearly  given  by  the  second 
time  derivative  of  the  incident  field's  pulse  amplitude. 

4.  TRANSPARENCY  CONSIDERATIONS 


The  transparency  of  a  single  carbon  or  graphite  fiber  In  an  epoxy 
matrix  Is  considered  based  upon  the  fiber  elastic  constants  reported  by 
Smith  (Reference  15).  Adopting  the  notation  in  Reference  15,  the  coordi¬ 
nate  axes  are  defined  with  the  3  axis  parallel  to  the  fiber.  For  a  plane 
dilatation  pulse  normally  Incident  upon  the  fiber,  e3  *  0.  Characteristics 
of  the  fiber  (modeled  as  having  hexagonal  crystal  symmetry)  imply  that 


’16  ** 26 


0,  c 


11 


"  c22*  and  c66 


’ll 


-  c12). 


Using  these  to  evaluate 


Equation  25  for  the  product  e^eje|ce]c^ »  the  amplitude  M^bgr  Is  ob¬ 
tained  according  to  Equation  34  as 
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Figure  8.  Comparison  of  an  Experimental  Result  (Solid  Line)  with  a 

Theoretical  Prediction  (Broken  Line)  for  Backscatter  from  a 
Cylinder  for  ka  ■  0.32 
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Note  that  this  Is  the  same  result  obtained  for  Isotropic  media  in  Equation 
39.  This  Is  a  consequence  of  the  hexagonal  model  of  the  fiber  and  the 
normal  Incidence  of  the  dilatation  pulse.  It  is  well  known  (Reference  16) 
that  a  hexagonal  structure  will  appear  the  same  to  a  normally  incident 
wave.  Independent  of  where  e  lies  in  the  1-2  plane,  l.e.,  the  fiber  is 
"transversely  isotropic." 

Table  1  illustrates  the  scattering  amplitudes  for  the  fibers  consi¬ 
dered  in  Reference  15. 

TABLE  1 

FIBER  CONSTANTS  AND  SCATTERING  AMPLITUDES 


Fiber 

P 

(g/cm3) 

C11 

(1010Pa) 

Mfiber 

(  "flberY 

(  "void  ) 

WYB 

1 .32 

3.13 

5.10 

6.50 

T-25 

1.38 

1.94 

2.89 

2.09 

T-40 

1.57 

1.42 

2.06 

1.06 

T-50 

1.67 

1.25 

1 .82 

0.83 

T-50S 

1.69 

1.22 

1 .78 

0.79 

T-75S 

1 .88 

0.97 

1.47 

0.54 

VYB 

1.53 

4.92 

8.69 

18.88 

PAN 

1.72 

2.31 

3.88 

3.76 

HTS 

1.69 

1.60 

2.51 

1.58 

T-400 

1.77 

2.33 

3.97 

3.94 

2 
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from  E 


p°  =  1.16  g/cm^  and 

3.27  x  10^  Pa  and  v°  =  0.35  according  to 


The  epoxy  density  was  taken  as 
o 


'11 


was  calculated 


*  EO  -  v) 

11  (1  +  v)  (1  -  2v) 


(121) 


The  last  column  In  the  table  compares  the  backscattered  energy  from  the 
fibrous  material  to  that  from  a  void  of  the  same  shape.  For  only  three 
materials  is  the  backscattered  energy  less  than  for  a  void. 


Practical  consideration  of  detecting  porosity  in  an  array  of  fibers 

must  also  include  the  geometry  of  the  scatterer,  i.e.,  the  "strength"  of 

the  scatterer  mentioned  in  the  section  on  long  wavelength  scattering. 

3  2 

The  factors  4irs  /3  and  irLc  are  what  is  to  be  considered ,  not  the  factors 
2  3  2 

4bLaQ  and  2itaQb  appearing  in  Equations  73,  95,  and  103,  which  are  due 
to  the  time  normalization  employed.  Consider,  for  example,  fibers  8 
microns  in  diameter,  illuminated  by  a  normally  incident  wave  which  is 
approximately  planar  over  a  length  of  1/4  inch  (the  size  of  a  typical 
transducer  face).  Without  considering  the  material  parameters,  the 
amplitude  response  from  a  single  sphere  will  exceed  that  of  a  single  fiber 
for  porosity  diameter  greater  than  about  3.3  mils  (1  mil  =  25.4  microns). 
By  Including  the  material -dependent  amplitudes  from  Table  1  and  requiring 
that  (4tts  /3)Mvojc|  be  greater  than  uLc  M^ber,  one  finds  that  the  lower 
limit  of  porosity  size  that  produces  a  larger  backscattered  response  than 
an  8  micron  diameter  fiber  varies  little  compared  to  the  range  of  energy 
ratios  in  the  table.  For  energy  ratios  varying  from  0.54  to  18.88,  the 
lower  limit  of  "detectable"  porosity  varies  from  3  to  5.5  mil  in  diameter. 
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SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 


1 .  CONCLUSIONS 

The  first  Born  approximations  to  solutions  of  a  time-domain  integral 
equation  were  used  to  obtain  the  backscattered  dilatation  wave  response 
from  spherical  and  cylindrical  inclusions  of  arbitrary  homogeneous  ani¬ 
sotropic  elastic  media  embedded  within  a  homogeneous  isotropic  host.  For 
relatively  large  values  of  normalized  center  frequency  ka,  the  time  wave¬ 
form  responses  from  cylinders  and  spheres  are  markedly  different  due  to 
the  difference  in  slopes  of  the  scatterers'  cross-sectional  areas  which 
the  wave  sees  as  it  first  meets  the  scatterer  and  then  as  it  exits  the 
scatterer.  For  ka  -*■  0  the  responses  from  the  cylinder  and  sphere  have 
identical  time  form  and  are  given  by  the  second  derivative  of  the  input 
pulse  time  profile;  they  differ  In  amplitude  and  the  difference  depends 
upon  the  volume  of  the  scatterer.  Excellent  agreement  with  an  experi¬ 
mental  result  for  scattering  from  a  cylindrical  void  was  obtained  at  a 
value  of  ka  -  0.32. 

A  "transparency  condition"  was  obtained,  which  states  that  for  certain 
combinations  of  both  density  and  stiffness  of  the  scatterer  and  host,  the 
scatterer  appears  transparent  to  the  incoming  wave  in  the  first  Born 
approximation.  For  inclusion-host  combinations  which  do  not  satisfy  the 
condition  exactly,  a  useful  quantity  is  the  energy  scattered  from  an 
inclusion  normalized  to  that  scattered  from  a  void.  This  is  of  practical 
significance  for  non-destructive  Inspection  of  fiber-reinforced  composite 
materials  with  long  wavelengths.  A  lower  limit  exists  on  the  size  of 
spherical  voids  which  give  a  larger  amplitude  response  than  a  cylindrical 
Inclusion  of  a  given  size  based  upon  both  the  relative  transparency  of 
the  Inclusion  and  the  volume  of  both  scatterers  illuminated  by  the  wave. 
For  an  8m  diameter  fiber  Illuminated  over  1/4  inch,  this  lower  limit  on 
the  spherical  void  size  ranges  from  3  to  5.5  mil  in  diameter  for  a  wide 
variety  of  fiber  compositions.  Since  these  values  are  within  the  range 
of  practical  Interest  (1  -  20  mil).  It  appears  from  this  simple  analysis, 
that  an  Inspection  technique  based  solely  upon  backscattered  waves  should 
not  be  dismissed  without  further  analysis. 
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2.  RECOMMENDATIONS 

Extensions  of  this  study.  Immediately  applicable  to  the  non-destructive 
Inspection  problem,  should  Include  considerations  of  backscatter  from 
spheres  In  the  proximity  of  an  Infinite  linear  array  of  parallel  cylinders. 
All  angles  of  incidence  with  respect  to  the  plane  normal  and  the  direction 
parallel  to  the  fibers  should  be  considered,  since  normal  incidence  alone 
Is  not  usually  used  in  current  inspection  techniques.  Investigation  of 
the  effects  of  periodic,  almost  periodic,  and  random  spacing  of  the  array 
could  prove  useful.  Since  the  number  of  cylinders  needed  to  simulate  the 
problem  numerically  would  be  very  large  in  the  long  wavelength  regime  of 
interest,  the  computer  program  developed  for  use  in  this  study  would  not 
be  useful,  providing  as  it  does,  complete  solutions  for  the  scattering 
from  each  fiber  in  the  array.  Further  analytical  work,  summing  the 
responses  from  the  cylinders,  should  be  carried  out  first,  and  the  scat¬ 
tering  from  spheres  obtained  against  a  "background"  response  from  the 
cylinders.  The  extension  to  many  randomly  distributed  spheres  near  the 
array  would  then  provide  a  model  which  closely  approximates  the  physical 
system  In  the  non-destructive  Inspection  problem. 

Basic  questions  regarding  the  validity  of  the  time-domain  first  Born 
approximation  need  to  be  answered.  The  statement  of  the  approximation. 
Equation  14,  suggests  that  perhaps,  larger  ka  problems  might  be  solved 
by  this  technique  provided  the  perturbations  in  density  and  stiffness 
presented  to  the  wave  are  sufficiently  small.  Certainly,  as  Ap  -*■  0  and 
Ac1jkl  **■  backscatter  goes  to  zero,  and  a  regime  of  first  Born 

validity  should  exist  for  small  material  perturbations.  Another  nagging 
point  Is  the  artificial  unboundedness  of  the  material -dependent  amplitude 
factor  M  In  the  ratios  of  density  and  stiffness.  An  upper  bound  needs  to 
be  obtained  for  the  magnitude  of  M  beyond  which  the  first  Born  approxima¬ 
tion  ceases  to  be  usable. 
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